ENERGY CONCENTRATION AND EXPLICIT 
SOMMERFELD RADIATION CONDITION FOR THE 
ELECTROMAGNETIC HELMHOLTZ EQUATION 



MIREN ZUBELDIA 



Abstract. We study the electromagnetic Helmholtz equation 
(V + ib(x)) 2 u(x) + n{x)u{x) = f(x), x G R d 

with the magnetic vector potential b(x) and n(x) a variable index of 
refraction that does not necessarily converge to a constant at infinity, 
but can have an angular dependency like n(x) — > rioo f]fjj as \%\ oo. 
We prove an explicit Sommerfeld radiation condition 
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for solutions obtained from the limiting absorption principle and we also 
give a new energy estimate 
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' , . dx < +oo, 
1 + \x\ 

which explains the main physical effect of the angular dependence of n 
at infinity and deduces that the energy concentrates in the directions 
given by the critical points of the potential. 



1. Introduction 

Let us consider the following Helmholtz equation 

(1.1) {V + ib(x)) 2 u{x) + n{x)u{x) = f(x), x£R d 

where b = {b\ , . . . , bj) : W 1 — > M. d is a magnetic potential and n : M. d — > R is 
a variable index of refraction that admits a radial limit 

(1.2) nix) —7- Tioo ( — ] as la^l — > oo. 

Km J 

In this work we are interested in the study of the existence and uniqueness 
of solution of the equation (jl.ip with an appropriate radiation condition 
using the limiting absorption method as well as some new estimates that 
characterize the behavior of the solution at infinity. 
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The radiation conditions are known to be necessary for the uniqueness 
of solutions to (jl.ip . as first remarked by Sommerfeld [So] in the free case 
b = = V . Moreover, there is a strong connection between this topic and 
the limiting absorption principle for Schodinger operators, to which a lot of 
research has been devoted (see for example, [El] . [15] . [Ku] . [A"] . [SI] . [Ml] . 
[M2] . [MU] . [S2], [I], [Bg], [E2], [5J, [PV2] . [RT] . 0). Let us now give a brief 
picture about the recent advances on Sommerfeld radiation. 

As a first result we mention Eidus [El], where it is showed that there 
exists a unique solution u(X, f) of the equation (jl.ip with n(x) = A + V(x) 
in M 3 satisfying the radiation condition 

2 



(1.3) 



lim 

r— too 



du 



d\x\ 



iX l > 2 



it 



da(r) = 0. 



Here bj(x) is assumed to vanish close to infinity and the electric potential 
satisfies V(x) = <3(|x|~ 2 ~ Q ) with a > g at infinity. In 1972, Ikebe and 
Saito [IS] extended the result by Eidus to electric potentials of the form 
V = Xp + Q, where p is long range and Q is short range, obtaining the 
precise radiation condition 

i\ 1/2 At u 



(1.4) 



(V + ib)u 



dx 



(i + N) 



l-<5 



< +oo, 



where < 5 < 1 is a fixed constant. This result has been recently improved 
by Zubeldia in [Z], by adding some singularities on the potentials at the 
origin and extending the range of 5 to < 5 < 2, in (jl.4|) . In the purely 
electric case 6 = 0, under similar assumptions on V, Saito proved in [SI] , 
[S2] another type of radiation condition, 



\Vu - i(VK)u\ 



dx 



< +oo, 



(i + \x\y- s 

being K = K(x, X) an exact or approximate solution of the eikonal equation 
(1.5) \VK\ 2 = A(l +p{x)). 

Here the principal order of VK has the form 

VK(x,X) = $(x,A)^. 

\x\ 

Notice that the unit normal to the sphere appears in the radiation con- 
dition, in all the above mentioned papers. 

In 1987, Saito considers more general long range potentials and proves 
the existence of a unique solution of the equation (jl.ip for n(x) = A(l + 
p(x)) + Q(x) with a nonspherical radiation condition 



(1.6) 



/ (V + ib)u - iVXVKu 

J\x\>l 



dx 



(i + M 



\X-S 



< +oo 



where VK is the outward normal of a surface which is not a sphere in 
general and satisfies the eikonal equation (|1.5p . More precisely, he consid- 



ers potentials p(x) = Xp(x) such that p{x) = 0(1), = 0(\x\ 1 ) and 
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q®.q x = 0(\x\ 2 ) at infinity. Under the same assumptions on the poten- 
tials, Barles establishes in [B] the existence of solution of the corresponding 
eikonal equation for \x\ > Rq, for Rq, A large enough. 

Some years later, in 2007, Perthame and Vega |PV2] showed that, in 
the purely electric case 6 = 0, if one puts some further restrictions on the 
potential n(x) as 

(x ■ Vn(x))- 



2 2_. sup ■ 



n(x) 



< 1, 



where C(j) denotes the annulus {2 J_1 < \x\ < 2 jf }, while (a)_ = — min{0, a} 
is the negative part of a G M and if there exists noo G C 3 (S d ~ l ) such that 
^cxd(w) > no > with 



(1.7) 



\n(x) - (a;) | < ^ooM- 



T > 0, n > 0, 



then the following precise radiation condition holds 

2 



V-u 



'W>i 

Note that the spherical term 
(1.8) 



n 



in 



1/2 



1/2 



dx 



< +oo. 



a; 
x 



appears in this formula instead of the gradient of the phase as in (|1.6j) . This 
apparent contradiction can be explained by the existence of some extra 
energy estimate already announced in [P V4] ) of the form 



(1.9) 

for some uj 



V tl ,n r 



1 + \x 



■dx < +oo, 



#t G S d 1 . Here we define 



d_ 



n{ijj 



\x\ V n 



and 



V _L n(x) = Vm(i) — —V r u(x), V r u(x) = d r u(x) := — • Vu(x). 



Estimate (H} says that the points where |V tt) n 0O (w)| vanishes on the sphere 
are the concentration directions for the energy |n| 2 . In other words, the 
energy is not dispersed in all directions but concentrated on those given by 
the critical points of the potential. Thus we see that in this case the behavior 
of the solution at infinity can be very different to the one exhibited by free 
solutions. 

The role played by the critical points of n^ was already pointed out by 
Herbst [He| . where is considered the case when n(x) = A + V{x) with 



V(x) 



< a < 2, Vx G M d \{0}. 
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This potential is also studied in |FG] , |GVV| and [HSJ , for the study of the 
counterexamples of Strichartz inequalities for Schrodinger equations with 
repulsive potentials and the existence and completeness of the wave operator. 

Let us introduce some notations. For / : R d -)• C, we define the norms 

\ 1/2 
\f(x)\ 2 ) 

\x\<R I 



\Ro 
and 




R 



1/2 / Nl/2 

+ U>/ m 2 

V J\x\<Rq J 



3~ 

where J is such that < R < 2 J and C(j) = {x G R d : V < \x\ < 2^ +l }. 
The norms |||/|||i and N\(f) are known as Agmon-Hormander norms. We 
drop the index Rq if Rq = 0, getting then the Morrey-Campanato norm and 
its dual in the sense that 



fg 



<\\\f\\\N(g). 



Moreover, we denote the magnetic gradient by 

V b := V + ib 
and its tangential component by 

|V^| 2 = |V 6 «| 2 - |VH 2 , V r b u = ^-V b u. 



We also introduce the magnetic field B{x) associated to the magnetic po- 
tential b(x) which is given by the d x d anti-symmetric matrix defined by 

B - m -w. ^=(St-H) JJ ' 1 

The tangential part of the magnetic field B is given by 

d 

B T := ^B(x), (B T )j = V ^B kj . 
\x\ ^ \x\ 

1 1 k=l 1 1 



This quantity was introduced by Fanelli and Vega [FV] and is related to 
singular magnetic potentials. See [FVj and [Zlj for more details. 



One of the main contributions of this work is to extend the new energy 
estimate to the magnetic case, inspired to Perthame and Vega [PV2J. Let 
us consider the magnetic Helmholtz equation 

(1.10) (V + ib(x)) 2 u(x) +n(x)u(x) + ieu(x) = f(x), e > 0. 



Then our goal is to show the estimate (|1.9p for the solution u of this equation 
using integration by parts. 

To this end, we first need to prove suitable a-priori estimates to the so- 
lution of the equation f 1 1 . 1 1) . On the one hand, one needs to control the 
Morrey-Campanato norm of the solution and its magnetic gradient. On the 
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other hand, we emphasize that the estimate for the tangential component 
of the magnetic gradient 

iVjhtf 



-dx < oo 



turns out to be fundamental. For this purpose, we will require that n{x) 
and the tangential component of the magnetic field satisfy the condition 

(1.11) 2E +2^1^ 

However, for the energy estimate it will be necessary to put some further 
restrictions to n(x) as in [PV2| and also to the magnetic field B. In fact, we 
assume that 

(1.12) ^2 2j \B jk \ 2 < oo. 
Moreover, it is required that 

(1.13) there exists ( ^- j G C°° (V" 1 ) , n M ( ^- j > n > 0, 



nix) — n c 

' 'x\ 



<noolAl^r, !■'>'). /? >U. 



and 

(1-14) , , 

\x\ J \x\ 

Note that from (|1.13p and (|1.14p it may be concluded that 

(1.15) | tt, | < C and n > for \x\ large enough. 

We may now state the first result of this paper. 

Theorem 1.1. For dimensions d > 3, we assume \1.11\ )- [1.14\ )- Then the 
solution of the Helmholtz equation A1.10\) satisfies, for R > 1 large enough 

(1.16) / Vu:n 00 (^-Y^<C(l + e)(N ( j 



\x \ \n 1 / 2 

for some constant C independent of e and n. 

This theorem is the natural generalization of the result by Perthame and 
Vega in [PV2J, which as far as we know was not known for first order pertur- 
bations of the Helmholtz equation. Nevertheless, it does not seem that these 
conditions on the potentials are sufficient to prove the limiting absorption 
principle for the equation (jl.ip . 

In order to get the limiting absorption principle for the electromagnetic 
Helmholtz equation with long range potentials that in particular include 
those which are homogeneous of degree zero, we will follow Saito [S]. Let us 
consider the equation 

(1.17) (V + ib{x)fu + n(x)u + Q(x)u = f, 

with n(x) = A(l + p{x)) where p, Q : M. d — > M. can be interpreted as elec- 
tric potentials. Under suitable assumptions on the potentials, we will prove 
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the existence of a unique solution of the equation (|1.17p satisfying a spe- 
cific Sommerfeld radiation condition together with some a-priori estimates 
of Agmon-H6rmander type. We will use multiplier techniques based on in- 
tegration by parts, inspired by [PV1| . [F] and [Zj. We work with potentials 
that decay as in [S] at infinity and the most important issue is that we allow 
singularities on the potentials at the origin. 

It is worth pointing out that the self-adjointness of the electromagnetic 
hamiltonian 

L = (V + ib(x)f + V{x) 

with V(x) = Xp(x) +Q(x) is necessary for the limiting absorption principle. 
For this purpose, we need to require some local integrability conditions on 
our potentials. In what follows we will always assume that 

(1.18) bj e L 2 oc , V e Lj oc , J V\u\ 2 < v J |Vu| 2 , < v < 1. 

As a consequence, it follows (see |Z1| . chapter 1 for more details) that L is 
self-adjoint in L 2 (IR rf ) with form domain 

D(L) = {/ € L 2 (R d ) : j \V b f\ 2 - j V\f\ 2 < oo}. 

We can now state the second main result of the paper. 

Theorem 1.2. Let d > 3, p G C 2 (M d \{0}), r > 1 and fi > 0. We assume 
that 

(1.19) \B jk (x)\ + \Q(x)\ < j-^, if \x\ > r„, 

(1.20) \Q(x)\ < if \x\ < r , < a < 2, 

(1.21) \d a p(x)\ < C*\x\~ a (\a\ < 2), 

for some c > 0, where a = (a±, . . . ,aa) is an arbitrary multi-index with 
nonnegative integers ctj (1 < j < d), \a\ = a% + ■ ■ ■ + a^, d a = d" 1 • • ■ d^ d 
and C* is a positive small constant (0 < C* < I). In addition, let c\ small 
enough and we consider 

(1.22) \B\<j^^ \x\<r , 0<a<2, 
if d = 3 and 

(1-23) \ B \<T^ \x\<r , 

\ x \ 

if d > 3. We also require that the magnetic potential satisfies the condition 

(1.24) |V • A\ < c\x\~ 2 , 

for some c > 0. Then, for any A € [Ao,Ai] with < Ao < Ai < oo, there 
exists a unique solution of the Helmholtz equation {l-llty satisfying 

(1.25) AIIHII^ + lllV^IH^^iVK/)) 2 
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and the radiation condition 
(1.26) 

/ V b u - iX^VKu 2 J X <C f (1 + \x\) 1+s \f\ 2 dx, 

J\x\>l + 6 J\x\>l/2 

for any < 5 < 1 such that 5 < fi, where C = C(Ao) and K is solution of 
the eikonal equation 

(1.27) \VK\ 2 = l + p(x). 

Remark 1.3. We point out that condition (|1.24j) is only needed for the unique 
continuation property, which is fundamental for proving the uniqueness re- 
sult (see Theorem 14.41 below). We use a unique continuation result proved 
by Regbaoui [R] . 

Theorem II .21 is the analog to the result by Saito in [S], generalized to pos- 
sibly singular potentials. Observe that in our approach it will be necessary 
to solve the eikonal equation (|1.27|) . 

Barles [B] proved that under the assumption (jl.21[) and for C* small 
enough, there exists a solution of the equation (|1.27p for \x\ > Rq with Rq 
large enough, see section [2 In general, one can not expect that the vector 
VK points at the direction A. An illustrative example given by Saito [S] is 
to consider 

p{x) = —\\A- 

Then p(x) satisfies (|1.2ip for A large enough and 

K{x) = a{\)\x\ -6(A)a?i 

with 

f a(A) = i[(l + l/A) 1 /2 + ( i_i/ A )i/2 ] 
\ & (A) = |[(1 + 1/A)V2_ (1-1/A)V2] 

is a solution of the eikonal equation (jl.27|) . This boundary condition differs 
from ours in all points except when Vn = (n = A + p). In the trivial case 
p(x) = 0, one can take K(x,X) = \x\. 

Note that assumptions needed to obtain the energy estimate and those 
for the limiting absorption principle are different and not comparable. On 
the one hand, if n = and regular, (jl.lip is trivially fulfilled. On the 
other hand, condition (jl.2ip with n = A(l +p) does not imply the existence 
of the limit Uoq. In addition, condition (|1.7j) does not need any regularity 
assumption on n(x) as in (jl.2ip . It is easy to see that if besides (jl.2ip . we 
require 

(1.28) \d r p{x)\ < c 2 \x\- 1 - f ", \x\ > 1, 

for some C2 > 0, \i > 0, then the index of refraction n{x) admits a radial 
limit rioo as \ x \ ~ * 00 • Moreover, it follows that 

\n(rui) — 7100(0;) | < r|x| _/i 
for r > 0, where r = |x| and u) = A. See [PV2| for more details. 



s 
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A combination of Theorem 11.11 with Theorem 11.21 will allow us to deduce 
an explicit Sommerfeld radiation condition. Let us consider the Helmholtz 
equation 



(1.29) 



(V + ib) 2 u + A(l +p)u = f 



and we denote n{x) = A(l +p(x)). The following result complements that 
of Saito [S] when Q = and extends the one given in |PV2j to the magnetic 
case. 

Theorem 1.4. For dimensions d > 3, assume U.ll\) and hl.21\) . Then for 
sufficiently small C* > and for any A € [Ao, Ai] with < Ao < Ai < oo ; 
there exists a unique solution of the Helmholtz equation H1.29\) satisfying 



(1.30) 



. 1/2/ \ 
VfeW — in 1 (X)-r-rU 

\x\ 



£<C S 
\x\ 



[ (l + \x\) 1+S \f\ 2 dx, 



for some 5 > 0. Moreover, if there exist noo, T > and fj, > such that 



n(x) - ( — 



< n(x)- 



for \x\ large enough, 



then it follows that 



(1.31) 



VfeM — in 



1/2. 



-u 



dx 



< C 



|/| 2 dx. 



Note that the spherical term 

\x\ 

appears in this formula instead of the phase as in (jl.26p . where VK is the 
outward normal of the surface |-fT(x,A)| = A, which is not necessarily a 
sphere. This apparent contradiction can be explained by the extra estimate 
(|1.16p on the energy decay which applies for the above example. In fact, 
it can be interpreted as a concentration of the energy along the directions 
given by the critical points of Uoq. In other words, the Sommerfeld condition 
hides the main physical effect arising for a variable n at infinity; the energy 
concentration on lines rather than dispersion in all directions. 

The rest of the paper is organized as follows. In the next section, we give a 
brief exposition of the eikonal equation and its properties that will be useful 
for the proofs of the main results. Section [3] will be concerned with the new 
energy estimate. We will prove Theorem 11.11 showing first the appropriate 
a-priori estimates given in Theorem 13. II that permits to deduce the desired 
conclusion. In section [3] we proceed with the study of the limiting absorption 
principle for the equation (|1.17p . We conclude the proof of Theorem 11.21 
following the ideas in [Z]. Section [5] provides a detailed proof of the new 
explicit Sommerfeld condition given in Theorem 11.41 The fundamental tool 
of the proofs are Morawetz-type abstract identities based on integration by 
parts, which are established in Appendix (see Lemmas 16.11 and 16. 3 [) . 
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Notation. Throughout the paper, C denotes an arbitrary positive constant 
and k stands for a small positive constant. In most of the C8.SCS, K, will come 
from the inequality ab < Ka 2 + j^b 2 , which is true for arbitrary n > 0. In 
the integrals where we do not specify the integration space we mean that we 
are integrating in the whole R d with respect to the Lebesgue measure dx, 
i.e. / = f Rd dx. 



2. The Eikonal Equation 

In order to determine the phase arising in the Sommerfeld radiation con- 
dition (jl.26p and to conclude the explicit one (|1.3ip . we need to solve the 
eikonal equation 

(2.1) \VK\ 2 = 1 +p(x), x£R d . 
Setting 

(2.2) g(x,C*) = \x\- 1 K(x,C*), 
we derive the following Hamilton-Jacobi equation 

(2.3) \g\ 2 + 2rd r g+\x\ 2 \Vg\ 2 = l+p(x), x£R d \{0}. 

From (|2.3p . under assumption (|1.2ip . Barles showed in [B] that there exists 
Co > such that for any C* < Co the differential equation f)2.1|) has a 
solution ip = ip(x,C*) for |x| > tq which satisfies 

(i) (p(x,C*) is a real-valued C 3 function for \x\ > vq. 

(ii) For any |x| > tq and C* < Co, 

(2.4) c < 5 (a;,C*)<ci 

with positive constants Co and c\. 

(iii) When C* -> 0, 

1, 3=0 
0, j = 1,2,3 



(2.5) \x\ J (&'g)(x,C*) — 

uniformly for x £ {x G R d : \x\ > ro}. 



Therefore, one can easily deduce the following identity that will be very 
useful in section HI 

Lemma 2.1. f[S], Lemma 2.5) For the solution K of the eikonal equation 
H2.1\) and for 1 < i,j < d, the following identity holds 



r» e \ Q2K \^ K ? , IdKdK 1 

(2 ' 6) 9^ = " l^'dx'i'dxj + K F ^ C ^ 

where Fij(x,C*) is a bounded function of x for \x\ > ro suc/i i/iai 

(2.7) lim sup \Fij(x, C*)| =0 = 1, . . . , d) 

|a!|>ro 

and 

1 i = j, 



i ^ j. 
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Proof. Setting K(x,C*) = \x\g(x,C*) and x^ = |n, we have 

dK , . dg 

— = XiQ + x — 

ax ? - oxj 



/ s dK ^ dg 

(2-8) =Xl5 + | x |^i. 



Then, 



■ i M ii x%Xj . . -|- x j t ij „ + \x 

oxiOXj \x\ \x\ oxj axj dxidxj 

which can be written as 

< 2 - 9 > ^-^-^ + ^^ 

with 



^L + x\x\— + \x\ 2 ^ g 
On the other hand, from (12.811 it follows that 



n ( n*\ — \ ~ I i i ~ i _i_ i 

<~r?.j \X)Lj ) — I Xj | X| — -+- Xj \ X\ — T" \X 



~ n — dK i i dg 
~~ dxi \ dj \dx l ' 

\VK\ 2 = g 2 + 2#|x|x • V# + |x| 2 |V#p 



Thus we obtain 



XiXj o 1 dK dK I f dg dK dK dg l2 9o <9<7 
-<? = — — — — — |x| — — — h |x| — — — |x| 



i<C X ctej dxj K \ dxi dxj dxi dxj dxi dx 3 

which together with (|2.9|) . gives (|2,6p with 

da dK 

(2.10) ^ = Gij - <%(2|x|x • Vg) + |x| 2 |V<7| 2 ) + N^^T 

, <91f <9o , , 9 dg dg 

+ x 1 ' 1 



9xj dxj dxi dxj 



The relation (|2.7p follows from ()2.5|) and the lemma is proved. 



□ 



In order to prove the explicit condition (|l,3ip we shall deduce (see section 
[5]) the estimate 

(2.11) / \V ± Ku\ 2 — Vt < +00. 

J l + |x| 

This is an energy estimate in itself which says that u concentrates along 
the critical points of VJ-ftT. In fact, from the hypotheses (|1.28p for the 
potential p(x), it follows that these critical point coincide with those of 
Voj^oo establishing a relation between the energy estimate (|1 . 16[) and (|2.1ip . 



Lemma 2.2. f[PV2], Theorem 3.2) Under assumptions il.21\) and fT. 

the solution to \2. 3\) satisfies for C* small enough and x ^ 

(2.12) \d r g\ < C*r- X -^, 

and g f r ]fy) ~~ ^ #00 fjfjj as r ^ 00, a smooth solution to the equation 

(2.13) g oa (u) 2 + \V u g oo (u)\ 2 =n 0B {u), u e S^ 1 . 
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Moreover, 

(2.14) \V^K\ = \V ugoo {u)\ + 0{r-») 
and 

(2.15) < CilV^I < iV^nool < C 2 |V w ^oo|. 

Observe however that for the energy estimate fll,16p we do not need the 
existence of a solution to the eikonal equation (|2.ip which could well not 
exist. 



3. The energy estimate. Proof of Theorem 11.11 

The purpose of this section is to extend the result by Perthame and Vega 
|PV2j to the magnetic case. To be more precise, we are interested in proving 
the energy estimate (|1.16p given in Theorem 11.11 for solutions u G H\{U. d ) 
of the magnetic Helmholtz equation 

(3.1) V 2 ti + n(x)u + ieu = f(x), e > 0. 

This estimate uses in a strong way the a-priori estimate for the Morrey- 
Campanato norm of the solution u of the equation (|3.1|) as well as the 
estimate for the tangential part of its magnetic gradient. 

Let us consider n(x) > such that 

(3.2) n = n\ +ri2 with ri2 G L°° , 

1/2 

(3.3) ||n 1 u\\ L 2 < (1 — co)||Vu||^2 for some Co > 0, 

, (x- Vn(x))_ + 2 2 i\B T \ 2 
3.4 2^sup^ l -^-:=P<l, 

n{x) 

where C(j) = {x G R d : 2 j ~ 1 < \x\ < 2 j } and (o)_ denotes the negative 
part of a G R. 

Then it follows the following result. 

Theorem 3.1. Let d > 3 and assume that i3.2\) - pT4\ ) hold. Then the 
solution to the Helmholtz equation H3. 1\) satisfies 

(3.5) M 2 :=|||V b n||| 2 + 111^/^1112+ nVhl] " 



< C(e+ \\n 2 \\ L oo) ( N 
where C is independet of e and n. 



nV2 



Proof. The proof is based on the identities which are stablished in Appendix 
and follows the same arguments of the proofs of Theorem 2.1 in [Z] or 
Theorem 1.1 in [PV2J. Thus we give only the main ideas of the proof. 

Let R > and we consider the functions ip and tp given by 



Vtp(x) 



^ if \x\ < R, 
iff ifN>«, 
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ip(x) 



2^ if \x\ < R, 
if \x\ > R. 



Let us add the identity (j6.6[) to (|6.2p with the above choices of the multi- 
pliers, respectively. Then, analysis similar to that in the proof of Theorem 
2.1 in [Z] gives 



(3.6) 



J\x\<R J\x\>R \ x \ ./|a;|=R 



< o / ($-")-M 2 + - / |x||B T ||V & u||u| + / |5 T ||V^u||u| 
+ e J \X7 b u\\u\ + 2 y |/||V 6 «| + C f 



The terms related to / can be treated as in the proof of Theorem 2.1 in 
[Z] . obtaining 



[ l ~lxT + [ l/HV6n|<J|||V fc n||| 2 + sup-^ f 



(3-7) / / |/||Vm/Uk( ii!V 6 «||| 2 + sup-^ / |// 

+ C K (iV(/)) 2 , 
where k denotes an arbitrary positive small constant. 

Let us study the potential terms. On the one hand, we have 



2 



(3 ,) ij^ U u?4eL 



(x ■ Vra)_ 



2J" 1 



n\u\ 2 



n 



~ ^ n 



On the other hand, let J such that 2 J 1 < R < 2 J . Then by Cauchy- 
Schwarz inequality yields 



R 







< 


si 








1 


< 






4R 



1 1 

\ 2 / \ 2 

2 1 I / I 1 2 1 o |2i„ ,|2 I 



(3.9) 

i 

/ |x| 2 |5 T | 2 |n 
J\x\<R 

\ I |VH 2 + E^^"'- 1/2 - 

K J\x\<R TT-, 



and 

1/2 / \ 1/2 



(3.10) f \B T \\v£u\\u\ < ( / 
J|x|>fl y7|a 



< 


t 




IvM 2 






'|x|>i? 




< 




/ 


IvM 2 




4 . 


l\x\>R 





b||B T | 2 |d 2 



+ / ? — \ — \\\ n u \\\ ■ 

j>J y ' 
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Finally, let us analyze the e term. In this case, the a-priori estimate (|6.5p 
reads as 

r \V b u\ 2 < [ n\u\ 2 + [ \f\\u\, 



which together with assumptions fj3.2[) - f|3.3[) implies 

J \y b u\ 2 <C^J n 2 \u\ 2 + J \f\\u\\ . 
Hence, by the same method as in [Z] it follows that 
(3.11) e j \V b u\\u\ < /€|||n 1 / 2 u ||| 2 + (7^(6 + sup |n 2 |) (n (-^\\ . 

As a consequence, plugging (|3.7p - (|3.1ip into f|3.6j) and taking the suprem- 
mum over R, by condition (|3.4j) we get (|3.5p . which is our claim. □ 



Remark 3.2. The dimension two is a special case. In this case, with the 
above choice of multipliers it follows that 

C 

A(2ip - Alp) < ~— pr. 

\ x \ 

Because of this singularity at zero, we cannot recover the full result (|3.5p for 
d = 2 and we cannot reach the right behavior close to 0. With some mod- 
ifications in the proof (see |PVlj . section 5 for more details) and assuming 
that n > riQ > 0, then in the two dimensional case it may be proved that 

— 1/2 

for Rq = n Q the solution satisfies 
IIIVHII 2 * + \\\n^u\\\% + f H>Ba ^ff < C(l + e) (N Ro ^ " * 

The homogeneity of the above estimate makes it compatible with the 
high frequencies (replace n by fi 2 n). Moreover, (|3.5p allows us to get the 
new energy estimate. As we have already said, the estimate of the tan- 
gential component of the magnetic gradient given in (|3.5p turns out to be 
fundamental. In order to get it, we need the smallness assumption given in 
(|3.4p . However, the condition (|3.4p is necessary and can not be relaxed to a 
Coulomb type of decay, even if smallness is added (see |PV2j . Appendix for 
more details). 

We may now state the main result of this section, which together with 
the above result proves Theorem ll.il Its interest relies on the bounds stated 
in Theorem 13.11 

Theorem 3.3. For dimensions d > 3, we assume Hl.ll\) - [1.14\ ) and use the 
notation of Theorem \3.1l Then the solution of the Helmholtz equation 13. 1\) 
satisfies, for R > 1 large enough 



x 



2 I |2 

TT<C 
\x\ 



(3.12) 

l\x\>R 

for some constant C independent of e and n 



M 2 + [N 



f 



n 



1/2 
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Proof. The proof consists in using the basic identity (|6.6p with a test function 
that depends on the behavior of n(x) at infinity. We choose R > 1 such that 
(|1.15p holds and define 



i/> q (x) = q 



n.~ 



x 



for some non-decreasing smooth function 



l(r) 



for r < 1 
r for r > 2. 



Let us put ip„ into (|6.6p . obtaining 



(3.13) 



1 



V b U ■ V b 1p q ■ V b U 



j,fc=l ' 



3? / /VV g • V b u - - / fAiP q u + e I Vip q - V b uu. 



We simplify the notation using q = q y^J , = n^ui). Observe that 
dijj q _ q'rioo x k , q dn c 



(3.14) 
and 
(3.15) 



dxk R \x\ \x\ duk 



A/ 9" ^^-1 9 A 

R z R \x\ \xr 



The left hand side of the estimate ()3.12p will come from the term 
\Vp-Vij) q \u\ 2 = j Vn-Vil> q \u\ 2 



I 



which can be written as follows 



(3.16) 



Vn • Vf/>Ju| 



doo 



2 I |2 

\u\ 



+ J drU—rioolul 

+ J q\x\V t (n — n c 
= h + h + h. 



doo 



The first term on the right-hand side of ()3.16[) is the one that gives the lower 
bound of what we want to control. By (|3.4[) and (|1.13p . we get 

/2>- C||n ^ lko ° |||n 1/2 n||| 2 . 
R 

On the other hand, after integration by parts, by the diamagnetic inequality 

|V|«|| < |V 6 w| 
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(see |LL| ) and by (|1.14p . we obtain 
h 



I w [n ~ n 



A I 1 2 i O ^ n oo, I, | 

ooi i A^noo 11 + I— |x|V|ti||u| 



C 
~R 





drioo 


2 H! 




duo 


\x 2 



+ 



R 



iv 6 «nr, 



for k > 0. 



Let us estimate now the remaining terms of the identity (|3.13p . A straight- 
forward computation gives 

V b U ■ D 2 ^ q ■ V b U = -p^ n oo|Vj,'li| 2 + TTTTl V ft"' u | 2n oo 



+ 23? 



Rx 



V7r drioo i 

r\x\ ixp j VbU a7 • Vb u 



\x\ z 



Thus since q', q" and — J^j are supported in the ball {|sc| < 2R}, 

by the Cauchy-Schwarz inequality it follows that the absolute value of the 
above terms in the corresponding integral are bounded by 

C|| n oo||c 2 



R 



|V>| 2 



+ ll|V 6 u|||' 



Moreover, by (|3.15p and (j!.15[) one can easily check that 

Re J V(A^g) • V b uu < CWriooWcs J + j^p ) l V & n IM 



< 



— |||n ' u 

R 



As far as the term containing the magnetic potential is concerned, first note 
that by (|3.14j) and the fact that 



yields 



d 



B T ■ V b u = B T ■ Vfu, 



dip q n x q'rioo , q dn. 



j,k=i 

Thus by (TTTB and (fLT2l . we get 

d 



R 



j,k=l 



j,k=l 



+ 



R { 

C||"-oo||ci 

R 



iv^V /2 



In 1 ' 2 !.! 



x\>R \ x \ 



n l l 2 u\ 



We now turn to analyze the terms containing /. On the one hand, by (j3. 15|) . 
we have 

/ 



l/l|A^IM<^/ a 



\x\>R 



M < p 



n 



1/2 



\n l l 2 u\ 
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On the other hand, from (|3.14p it follows that 

|/||ViMV 6 u| < I ^|/||noo||V 6 u| + J ^||/| 



du 



C firm fl . T 
< 11 °° IK -' AT 

R 



n 



1/2 



IIIVHII 



\x\>R \ x \ 



Finally, the last term to be bounded is 

e j VV 9 • V b uu < C||n g llcl £ /" |V 6 u|M, 

which can be done as in ()3.11|) . 

Therefore, from the above inequalities, taking k small enough yields 



I 



du 



2 N! < ^i 

\x\ 2 - R 



Infill 2 + |||V b n||| 2 + 



|V^| 2 



\x\>R PI 



+ R \ UV2 



which gives (|3.12p and the proof of the theorem is over. 



□ 



A combination of the above two results asserts the desired energy estimate 
(|1.16p and the proof of Theorem 11,11 is completed. 

Remark 3.4. From Remark l3.2l the same result holds for the two dimensional 
case. 



Remark 3.5. Condition (jl . 14)) can be largely relaxed if, for example, n — n c 
is radial. It can be instead assumed the alternative conditions 



(3.18) 



n(x) — n c 



< n-, 



for \x\ > Rq, r > 0, 5 > 0, Rq > 1 



and that there exist $ < 1, S > and f > such that 



x\V ± (n - iioo) 



du 



dn r 



du 



+ n(x) 



ki- 



ln particular, when n — is radial then (|3.18p is sufficient. 

Remark 3.6. Note that in order to prove the energy estimate we impose 
conditions in each component of the magnetic field Bjk and not in the tan- 
gential component of B, as in the first result. This is due to the fact that 
the test function chosen in the proof of Theorem II .11 is not radial (see (|3.14p 
and (|3.17p above). 
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4. Limiting absorption principle. Proof of Theorem 11.21 

Our next goal is to prove Theorem 1 1 . 2 1 which asserts the limiting absorp- 
tion principle for the equation (|1.17p . following [S] and [Zj. In addition, the 
result will be true for short range electric potentials that can have singular- 
ities at the origin and more importantly, critical singularities at the origin 
for the magnetic field can be considered. 

To do this, let us consider the electromagnetic Helmholtz equation 

(4.1) Vlu + A(l + p)u + Qu + ieu = f. 

We first prove the corresponding Sommerfeld radiation condition and a- 
priori estimates for the solution u G H\(W i ) of this equation for AG [Ao, Ai] 
with < Ao < Ai < oo and e > 0. We next turn to show the uniqueness 
result related to this equation. Indeed, we will see that if u satisfies (|4.1|) 
with e = and / = 0, then u = 0. Consequently, we will be in a position 
to construct the unique solution of the equation (|1.17p with the radiation 
condition (|1.26|) at infinity as the limit of the solution of the equation (|4.ip 
when e — > in some sense. The detailed proof of this construction is given 
in [Z], see subsection 2.4. Thus we will omit it. 

Since the proofs are adapted from the ones of the main results of [ZJ, we 
will mainly focus on the analysis of the new terms, that is to say, p. 



4.1. Sommerfeld radiation condition. We begin by proving the Som- 
merfeld condition in terms of the Agmon-Hormander norm of the solution. 
This result may be proved in much the same way as Proposition 2.6 of [ZJ. 

Proposition 4.1. For dimensions d > 3, let Ao > 0, e > 0, / G L\ +6 and 

r-r-r-*, ~2 ~ 

assume that M.19\) holds. Let K be a solution of the eikonal equation \2.1\) . 
Then, there exists a positive constant C = C(Ao) such that for A > Ao and 
C* small enough, any solution u G H^iW 1 ) of the equation J^.ip satisfies 
for all R\ > ro 

(4.2) / \V b u - iV\VKu\ 2 ( 1 + 6(1 + K) s 

VK\ 2 \V b u\ 2 - \VK- V b u\ 2 
(1 + K) 1 - 5 

<C(l + e )f||H||? + (iV 1 (/)) 2 + / (1 + K) 1+S \f\ 2 



+ (1-6) 



Proof. The proof will be divided into three steps and it consists in the con- 
struction of the Sommerfeld terms which contain the square 

(4.3) \V b u - iV\VKu\ 2 = \V b u\ 2 + A|VA"| 2 |u| 2 - 2^V\VK ■ V b uu. 

We use the identities proved in Lemma 16.11 and Lemma 16.31 In this case, 
one must choose the multipliers depending on the solution of the eikonal 
equation. By abuse of notation, we write instead of a vector field E. 



9(r) 
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Let R x > r . We take a cut off function G C°°(R) such that < 9 < 1, 
d0/dr > with 

1 if r > jRi + 1 
ifr<i?i, 

and set 9{K) = 9(K(x, C*)). We define * : R -»• R such that 

^'(r) = (1 + r) 5 , 0<5<1 
and we set 9(K) = V(K(x,C*)). 

Step 1. Our first goal is to obtain the term 

J \VK\\\V b uf + \\VK\M 2 ) { / + { Ky-6 - 
Let us first compute 



flEEj + (E2J), 
with the following choice of the multipliers 

E = Vtp = y'(K)VK9(K) 

respectively. Let us analyze all the terms of the resulting identity by the 
same method as in the proof of Proposition 2.6 in [Zj. In what follows, k 
denotes an arbitrary positive small constant and we use the same letter C 
for any positive constant. 

On the one hand, by (|2.6p and the fact that 9' is nonnegative, we get 

+ / 0(K) (&±*l - ^p,) {\VK\ 2 \V b u\ 2 - \VK • V b u\ 2 } 
+ / (1 E (V6) fc uF fcj (V^(iif) 



= h + I 2 + h, 

where 

h>{l-5)j {1 9 + { ^ 5 {\^K\ 2 \V b u\ 2 - \VK • V b u\ 2 }. 

On the other hand, observe that in order to get the term related to \u\ 2 
of the Sommerfeld square \V b u - i\ 1 / 2 VKu\ 2 , we need to use the eikonal 
equation (|2.ip . Indeed, we have 

. |2 5 [ \VK\ 2 \\VK\ 2 \u\ 2 „. T ^ 

ip X(l+p)\n\ 2 = -J { l + l K)l l J 0(K). 

Moreover, by the eikonal equation p(x) = \VK\ 2 — 1 and (|2.6p . it follows 
that 

("> £ = 2 E>^ *«-i* = i,...,* 

fe i=i J 
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Thus the other term involving the potential p gives 



k,j=l 

Let us treat now the terms containing the magnetic field B and the potential 
Q. Since c < |V-fT| 2 < c for some c, c > 0, by (jl . 19|) and Cauchy-Schwarz 
inequality we get 

S / §^ k B kmu(V b ) m u < C J \B km \\V b u\\u\(l + K) 5 e(K) 



k,m=l 



<kJ \VK\ 2 \V b u-i^/XVKu\ 2 -^- 



9(K) 



+ C K (v / A + l)||M||f. 
Similarly, by (j!.19[) we have 

|Vif| 2 |V 6 u - iyf\VKu\ 2 9{K) 



(1 + K) 



1-5 



Nil?. 



+ C7 K (VA + 1) 
In addition, since 

= y"{K)\VK\ 2 9{K) + y'(K)AK6(K) + #'(.K')|V.K'| 2 0'(.K') ) 
by (J23I), dS5D it may be concluded that 



- J vQ\u\ 2 + l - J QAip\u\ 2 < C||H||f. 
As a consequence, we get the inequality 
(4.5) S - J |VK| 2 (|V,n| 2 + A|V^| 2 |n| 2 ) (i ^ ) ) 1 _, 

+ (1-^)/ {l l { ^ s {WK\ 2 \V b u\ 2 -\VK.V b u\ 2 } 

- elm j 6(K)V'(K)VK ■ V b uu < -I 3 + h 

n f \VK\ 2 \V b u-iV\VKu\ 2 9(K) r- 2 

+ 2k J J \ 1+K) i-s ; +c(v^+i)||H|| 2 

-to I f U>'{K)VK- Vyu+ ^'{K)Ak) 6{K)u 

- | J y'(K)\VK\ 2 e'(K)u. 

Step 2. In order to obtain the desired square, the next step is to get the 
cross term 

9(K) 



-29 J V\\VK\ 2 VK- V b uu— 



+ K) 



IS' 
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Let us add to the above inequality (j4.5j) the identity f|6.3j) with the choice 
of a test function 

<p{x) = V\\VK\ 2 {1 + K) S 6(K). 

Hence, it follows that 

3 J Vp-V b uu = 5%V\ J +Ky-8 1 V-g| 2 ViT • V b uu 

+ J \VK\ 2 6'(K)(1 + K) S VK -V b uu 

J 1. 2 1 .7 



= h + h + h- 

The term I5 is used to complete the square can be 

upper bounded by 



k j \VK\ 2 \V b u - iV\VKu\ 2 ^ ^*y_ s + C K>Rl {l + 



A)||H 



In addition, denoting (Dx)iU = (V b )iU — iV~^§§-u, by (|2.10|) it may be 
concluded that 



-h + h + I 7 



< CC* 
Therefore, we deduce 



/ (1 + jf )S E {D K ) k uF kj (D K ) jU e{K) 
J k,j=i 

VK\ 2 \V b u - i^f\VKu\ 2 6(K) 



(l + K) l - & 

0{K) 



S - J \VK\ 2 \V b u- i V\VKu\ 2 {i + K)1 _ s 
+ {1 d + { Ky-s i\ VK \ 2 \ V "u\ 2 -\VK -V b n\ 2 } 

+ e^/\ j \VK\ 2 (1 + K) s \u\ 2 6(K) - e% J 9(K)(1 + K) 5 VK ■ V A uu 

<c(A + i)iirf + ( 3„ + «n 

-sft J f(l + K) s VK ■ (V^u + iX 1/2 VKu)9(K) 

-~ J fV'(K)(AK9(K) + \VK\ 2 9'{K))u. 

Step 3. In order to conclude the proof of Proposition I4.ll we build the 
Sommerfeld square (j4.3j) for the e term. 

Let us subtract the identity (|6.'2|) multiplied by e to the above inequality 
choosing the test function 
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so that we get 

^2 J \VK\ 2 (1 + K) 5 \V b u- i\ 1/2 VKu\ 2 6(K). 

In order to complete the estimate, by integration by parts and the a-priori 
estimate (|6.4|) . we have 



eft J V^p ■ V b uu = ~ J A(p\u\ 2 

<Ce f \u\ 2 < CWi(/)|||u| 
Furthermore, by (|1.19|) we deduce 

eJ,QH 2 <§- 2 J ' ^ 



\u\ 2 



\>r ^+\X\) 1 +H- S 

<Ce\\\u\\\l 

Finally, let us estimate the terms containing /. On the one hand, we have 



-ft J f(l + K) S VK ■ (V b u + i\ 1/2 V Ku)6(K) 



By (J23j), we get 

ft 



(l + K) 1 - 6 
+ C( K ) J(l + K) 1+s \f\ 2 6(K). 

J ty'(K){AK6(K) + \VK\ 2 6'{K))fu 

<c(llMlli + j(i + K) 1+5 \f\ 2 e(K)y 



By the a-priori estimate (|6.4p . yields 
M J (l + K) s fu6(K) 



<(x/i/i 2 (i+^) 1+5 w) 1/2 ( £ /h 2 ) 1/2 

<c(eJ |/| 2 (1 + K) 1+8 6(K) + HMD? + (7Vi(/)) 2 ) . 

Consequently, taking k > and C* small enough, we obtain (|4.2p and the 
proof is complete. □ 

Corollary 4.2. Under the assumption of Proposition \4^\ the solution u G 
H\(R d ) of the Helmholtz equation fr4.1\ ) satisfies 

\V h u-i\ l,2 VKu\ 2 

\x\>ro (1 + \ x \)^ J\x\>ro 



(4.6) / ' V+jdU-r +el (l + \x\) S \V b u- i ^ 2 VKu\ 2 



<C(l + e) \\\u\\\i + (N 1 (f)y + 



[ (1 + M) 1+, |/| 2 1 

J\x\>r J 
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for A > Aq, C* small enough and C = C(Xq). 



Proof. We need only take R\ = c^r® with co,ro given in section [2] and use 

(El). □ 



4.2. A priori estimates for A G [Ao, Ai]. Using the previous result, we are 
now in a position to prove the a-priori estimates for the frequency A varying 
in a compact set. We will deduce them by a compactness argument already 
used in [5] and [Z]. 



Proposition 4.3. For d > 3, under the hypotheses of Proposition \4-l\ let 
Ao > 0, A G [Ao,Ai], with Ai > Ao and e G (0, ei). Then, the solution 
u G H\(R d ) of the Helmholtz equation satisfies 

(4-7) A|||«|||f + |||V t «|||? < <7(l + B )(JVi(/)) a , 

where C = C{\q,E\). 

Proof. The proof is a combination of the proof of Proposition 3.1 in [SJ and 
the proof of Proposition 2.10 in |Z]. 

Let Bt be the interior of the closed surface Ey = {x : K(x,C*) = T} 
with T > ro and C* < Co, where ro and Co are given constants related to 
the assumptions of the potentials and the solution to the eikonal equation, 
respectively. Let us multiply the equation (|4,ip by u, integrate over Bt and 
take the imaginary part, obtaining 

vk _ r , , 2 



V b uu + e / \u\ =9 / fu. 



From this it follows that 



(4.8) 2VA9 / t^|- • Vb-uu < 2^9 / fu. 

iEj |viv| j Bt 

Let us integrate now the identity 

+ A|VK||-u| 2 = — !— \V b u - i^/\VKu\ 2 + 29VA 7 ^ r • V b uu 
\vK\ \vK\ Viv 



over the surface Et- Then by (|4.8p we get 

(4 - 9) 1 ( w + A|Vif ""' 2 ) £ L wk\ w " u ~ ! ^ Wu|2 

+ 2VLvi(/)||H||i. 

Let R > where p > ro, being co, c\ as in ()2.4p . Let us multiply both 
sides of (|4,9p by and integrate from pco to i?ci with respect to T. Hence, 
as |V-KT| 2 is lower bounded by a positive constant we have 

^[ (X\u\ 2 + \V b u\ 2 ) < ^ f \V b u-i\ l ' 2 VKu\ 2 

R Jpco<K<Rci J pcQ<K<Rci 

u\\\. 
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On the other hand, observe that since K = \x\g and cq < g < ci, yields 

{p < \x\ <R}C {pc <K< R Cl } C (— < \x\ < —} . 

[ ci c J 



Consequently, denoting j and j x by 2 jo ~ 1 < ^ < 2- 70 and 2 jW < < 
2 n , respectively, we deduce 

n 



4/ (A|n| 2 + |V 6 n| 2 ) < [ \V b u-i\ l ' 2 VKu\ 2 

RJp<\x\<R K ~" jo Jc(j) 

+ K\\\\u\\\l + C{n){N l {f)) 2 . 



Now, note that we are in the same position as in (2.42) of the proof of 
Proposition 2.10, [Z\. Therefore, by (jl .26|) . repeating the same reasoning to 
this case, it may be concluded that for R > 1 

1 / (A|^| 2 + |V 6 ^| 2 ) < ^||| u |||2 + C '(l + e)(iV 1 (/)) 2 . 

K J\x\<R 2 

Thus taking the supremum over R, the proposition follows. □ 

4.3. Uniqueness result. This paragraph deals with the uniqueness of so- 
lution of the equation (|1.17p . Let us consider the homogeneous Helmholtz 
equation 

(4.10) Vlu + X(l+p)u + Qu = 0. 

Then we formulate the uniqueness theorem as follows. 

Theorem 4.4. Let d > 3, Ao > and assume hi. 19]) . {1.2$ . Let u be a 

solution of the equation f^.iO ) with u,X7bU £ Lf oc such that 



(4.11) liminf/ (\V b u\ 2 + X\u\ 2 )da(x) -> 0, 

J \x\=r 



as r — f oo, 

'1*1= 



for A > Ao- Then u = 0. 

Moreover, if for some 5 > the condition 

(412) i,J Vi!, - iAl/2vA "" |2 (TW 

is satisfied, then \Jj-ll ) holds. 



1 < 00 



Proof. The proof follows by the same method as in the proof of Theorem 1.5 
in (ZJ. Although the analysis of the terms related to p are slightly different, 
the same conclusion can be drawn for this case. We give only the main 
ideas. For a fuller treatment we refer the reader to [Z] . 

The first step is to show 



Id 2 . 



(4-13) / (|VH 2 + |u| 2 )<^ / 

For this purpose, we multiply the equation (|4.10|) by the combination of the 
symmetric and the antisymmetric multipliers Vi^ ■Vf ) u+ ^Atpu + ipu and we 
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integrate it on the ball {\x\ < R\} for some R\ > R > r$. Then, we define 
a cut off function 6 with 



6(r) 



1 if r > 1 
if r < I 



and 8' > for all r. Set 9r[x) = 6 (JrJ an d for R > > 1 choose the 
multipliers as follows 

W(*) = ~0 R (x) 

and 

We do the computations for R large enough and we pass to the limit in R\. 
The next goal is to prove that for R > 2tq > 1 and any m > 0, then 



f \x\ m (\V b u\ 2 + u\ 2 ) < +oo. 

J\x\>R 



\x\>R 

This can be easily shown by induction. 

Our last claim is to prove the exponential decay of the solution. We define 
the test functions 

VV(x) = \x\ m+1 ^9 R (x), 
\x\ 

ip{x) = ±\x\ m R (x), 

for R > 2tq > 1 and we put them into the identities (|6.6|) and (|6.2|) . respec- 
tively. We add both equalities and analysis similar of that in [Z] , shows that 
for R large enough and for any t > 1, < 6 < 2/3, A > Aq, it follows that 



u\ 2 <C s e- tRS [1 + X + 



\x\>2R 



~R 



being C$ independent of t. Thus letting t — > oo, we obtain that u = almost 
everywhere in {\x\ > 2R}. The unique continuation property ([Rj) implies 
then u = almost everywhere in M. d . 

In order to deduce (|4.1ip from (|4.12p , first observe that solutions of (|4,10p 
satisfy 

• Vfeim = 0, 

-'T 



just multiplying the equation by u and integrating over Bt, the inside of 
the closed surface = {x : K(x, C*) = T}. Hence, we have 

f (\V b u\ 2 + X\VK\ 2 \u\ 2 )da(x) = f \V b u - iVXVKu\ 2 da(x), 
which together with (|4.12p gives (|4.11|) . 

□ 
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5. Explicit radiation condition. Proof of Theorem 11.41 

This section establishes the relation between the energy estimate (|X . 16[) 
and the Sommerfeld condition (|1.26p . We will see that when the variable 
index of refraction has the form n{x) = A(l + p(x)) and an angular depen- 
dency like n(x) — > noo (jjj^) as \ x \ ~~ > 00 > then the Sommerfeld condition 
(|1.26p at infinity still holds under the explicit form (|1.3ip . 

Proof of Theorem \l-4\ The proof follows |PV2j . Let us first recall the tan- 
gential estimate 

(5.1) / < C{N{f)f 



\x\ 

proved in Theorem 13.11 above and observe that (jl.26p provides 

(5.2) J \V b u - i\ 1/2 VKu\ 2 -^-^ <C J(l + \x\) 1+5 \f\ 2 dx. 

Hence, just looking at the tangential part of the above inequality, by (|5.ip 
it follows easily that 

(5.3) / A|V^n| 2 -^ < C [ + C A 1 + M) 1+ *l/| 2 

J\x\>r i- + \X\ J 1 -+- \X\ J 

<C j\l + \x\) 1+s \f 2 \. 

Furthermore, since n = A(l +p), from the eikonal equation ()2.ip we have 
n-\\d r K\ 2 = \X 1 / 2 V ± K\ 2 . 



Now, according to the properties (|2.5p related to VK, it is easy to see that 
d r K = g{x) + 0(C*) > 0. Thus we obtain 

(5.4) |,V^-„V^_W^< C | A V W 

In addition, looking at the radial part in (|5.2p we have 

(5.5) / \V r b u - iX^drKul 2 —^ <C f(l+ \x\) 1+S \f\ 2 . 

J\x\>r 1 + \ x \ J 

Consequently, by (|5.3p . (|5.4p . (|5.5p and the fact that 



\x\ 



x 



< \V T b u - iV\d r Ku\ 2 + \V\d r Ku - n^u\ 2 + |V^u| 2 , 



we get (|1.3U|) which is our first claim. 

Finally, assuming \n — I < C(l + \x\)- 5 and using (|3.5|) we conclude 



that 

2 



/ V.u-in^^u —l—<c[\f\ 2 (l + \x\ 
J \x\ l + \x\ J 

and the proof is complete. 



□ 
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6. Appendix 

In this section we state the key equalities that have been used in the 
proofs of the main results of this work. 

These integral identities are obtained by the standard technique of Morawetz 
multipliers, using integration by parts (see [FJ, Lemma 2.1. and |PV1| . 
Lemma 2.1). In order to carry out the integration by parts argument below, 
we need some regularity in the solution u. In general, it is enough to know 
that u G H\(R d ). See [Z], Appendix for more details. 

Let us consider the electromagnetic Helmholtz equation 

(6.1) (V + ibfu + A(l +p)u + Qu + ieu = f, A,e>0. 

Lemma 6.1. Let ip : R d — > R be regular enough. Then, the solution u G 
H\(R d ) of the Helmholtz equation \6. 1\) satisfies 

(6.2) 

J ip\\u\ 2 - J <p\V b u\ 2 + j (p(p + Q)\u\ 2 - K J Vp-V b uu = ^ J <pfu, 

(6.3) e j p\u\ 2 — 9 j Vip-V b uu = Q J ffu. 

Remark 6.2. Note that if we take ip = 1, then we obtain the following a-priori 
estimates 

(6.4) e j \u\ 2 < j \f\\u\ 

(6.5) J \V b u\ 2 < J(\ + p + Q)\u\ 2 + J \f\\\u\, 
that have been very useful throughout the paper. 

Lemma 6.3. Let ip : M. d i — > M be regular enough. Then, any solution 
u G H\(R d ) of the equation \6.1\) satisfies 

(6.6) j V b u ■ D 2 ip ■ V> + V(A^) • V b uu + J Vip ■ X^uu 
-3 / p^ B kj(Vb)jUu-\ [ AiPQ\u\ 2 -^ f QVifj-V b uu 

where D 2 ip denotes the Hessian ofip. 
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